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Рассмотрена задача моделирования гравитационного потенциала сосредоточенной массы, 
заглубленной в многослойной литосфере. 
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Введение 

Задача определения природы источника, вносящего аномальный вклад в 
распределении отклонения значений гравитационного потенциала от регионального 
распределения Ух, у, <), является некорректной, в том смысле, что измеренное реаль- 
ное гравитационное поле полевым гравиметром в данной малой области может быть 
обеспечено в общем случае значительным числом вариантов распределения источника 
аномального поля под данным гравиметром. Вместе с тем, практически задача в таком 
случае обычно ставиться в форме поиска «некоторого эквивалентного» геометрически 
простого тела (в нескольких вариантах), для которого ожидаемое трехмерное 
распределение гравитационной аномалии имеет одну из базовых модельных форм. 

Постановка проблемы 

В данной работе предлагается достаточно практичная модель расчета распре- 
деления масс исследуемой гравитационной аномалии, сводящаяся к конечномерной 
модели «элементарных» кубических объемов на трехмерной сетке, для которых можно 
указать несколько наиболее вероятных прогнозных распределений аномалий масс. 
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Анализ последних исследований 

Из практики грави-геологической разведки слабообследованных территорий в 
данной гравитационной разведке дополняется результатами приборного обследования 
данной аномальной зоны дополнительными технологиями (магнитовариационная 
разведка, геохимическая разведка пробы воздуха и водной компоненты грунта в зоне 
аномалии). Данная статья стимулирована практической задачей отыскания упавших 
метеоритов или признаков вулканической деятельности под слоем льда, несущих 
геофизические и геохимические эффекты как в толще снега или льда, так и вблизи 
поверхности ледового покрова однородной среды (типа снег и/или лед, характерные 
для Антарктиды). 

Данные исследования проводятся в результате замысла авторов организовать и 
выполнить экспедиционные работы в Восточной Антарктиде с намерением получить, 
по возможности, непосредственную регистрацию: 

[) ожидаемых признаков корреляции (в пространстве во времени) активизации 
тектонической и вулканической деятельности на Антарктическом континенте; 

2) наблюдаемых признаков ускоренного таяния некоторых ледников на 
локальных территориях континентального шельфа. 

По результатам экспедиции предполагается обеспечить информационную базу 
для получения средств прогнозирования зон и периодов ожидаемой активизации 
тектонической и вулканической деятельности на ряде территорий Антарктического 
континента. Планируется использовать для этих целей признаки возможной общности 
генезиса «вулканических популяций» Южной Америки и Западной Антарктики, на что 
указывают, в частности, недавние результаты глубоководного бурения океанического 
дна пролива Дрейка, а также ожидаемая прогнозируемость слабых землетрясений 
земных пород, покрытых антарктическими ледниками. 

Авторы [1] излагают свои доводы о целесообразности использовать в 
планируемых экспедиционных работах, в частности, сейсмографы и новейшие 
технологии подводных исследований, разработанные в Институте кибернетики имени 
В.М. Глушкова НАН Украины (группа ученых под руководством д.ф.-.м.н., проф. 
Писаренко В.Г., отдел математических проблем прикладной информатики - № 265) с 
применением оригинальных подводных автономных аппаратов Консорциума 
«ВаКоБойс$», имеющих видеосвязь и специальные технологии искусственного интеллекта [1]. 

Цель исследования 

При наличии данных о гравитационном потенциале, на поверхности Земли можно 
найти аномальные зоны, под которыми находятся области с иным литосферным составом. 

Основной материал 

В задачах гравиметрии часто используется скалярная величина: 


ат 
У=б| ==. 
У 
Алгоритм расчета поля гравитационного потенциала использует интеграл вида [2]: 


У-б] ат 
(хх О о ал 


(р) 
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Расчетная модель представляет собой трехслойный параллелепипед, площадь 
основания которого составляет 20х20 кв. единиц. Высота каждого слоя задается в 
условных единицах. Пользователь имеет возможность выбрать состав слоя, тем самым 
изменив соответственно его плотность. 

Алгоритм. Потенциальное поле рассчитывается для двумерной сетки {Х, У}. В 
каждом узле вычисляется трехкратный интеграл вида (1), где положено <, =0, с 


помощью трех вложенных циклов. Первые два цикла (по 1, ]) соответствуют подсчету 
потенциала одного глубинного слоя. Внутренний цикл по К рассчитывает вклад в 
потенциал от единичного вертикального столбца в узле двумерной сетки. В 
зависимости от высоты меняется значение плотности (числитель в (1)). Если в узле 
трехмерной решетки {1, ], К} находится аномалия, то нужно задать соответствующую 
плотность. После вычисления значений потенциала строится график трехмерной 
поверхности потенциала (рисунки 1-7). 

Программное обеспечение разработано в системе МАТГАВ [3]. 

Аномалия моделируется с параллелепипеда с поперечным сечением в | кв. ед. 
протяженностью 5 ед., залегающего на заданной глубине в центре расчетной области. 
Состав аномалии представляет одну из 6 железных руд (таблица 1). 


Таблица 1. Плотности железных руд 


№ Название руды Плотность, т/м? 
1 Магнетит 4,90 
2 Гематит 5,26 
Э Гидрогематит 4,28 
4 Гетит 4,20 
5 Лимонит 3,30 
6 Сидерит 3,96 
`& Гравитационный эксперимент ЕЛЕ Е 
Ее Ед \Мем пзеё Тоо5 ОебЖюор  МИпдом Нар м 


сен к <зоэги- а паес 


Параметры 


Без аномалии | С аномалией [ Разница| 
Глубина 1-го слоя 5 


Глубина 2-го слоя 10 
Глубина 3-го слоя 5 
Состав 1-го слоя Лед у 700 
Состав 2-го слоя Лед а 
Состав 3-го слоя Лед ый 600 
Глубина аномалии 15 


500 


Потенциал 


Состав аномалии |магнетит х 


400 
20 


Ширина оо Длина 


Рис. 1. Гравитационное поле над толщей льда. 
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Рис. 2. Изменение гравитационного поля над толщей льда, вносимое аномалией на 


глубине 15 ед. 
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Рис. 3. Изменение гравитационного поля над толщей льда, вносимое аномалией на 
глубине 10 ед. 
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Рис. 4. Изменение гравитационного поля над толщей льда, вносимое аномалией на 
глубине 5 ед. 
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Рис. 5. Изменение гравитационного поля над толщей базальта, вносимое аномалией на 


глубине 5 ед. 
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Рис. 6. Изменение гравитационного поля над толщей базальта, вносимое аномалией на глубине 10 ед. 
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Рис.7. Изменение травитационного поля над толщей базальта, вносимое аномалией наглубине 15 ед. 
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Выводы 

В ходе вычислительных экспериментов было выявлено, что отклонения 
гравитационного поля зависят от глубины аномалии и соотношения ее плотности с 
плотностью слоя, в котором она находится. 

В качестве примера был произведен расчет «надежности распознавания» 
метеорита плотностью 5 т/мЗ (плотность метеоритов колеблется от 3,5 до 8 т/мз) 
приближенной кубической формы размером в одну кубическую расчетную ячейку. 
Производился контрольный расчет обнаруживаемости предмета с аномальной 
плотностью 1,6 от плотности среды и (2, 3, 4 слои) с аномалией массы в 2.35-М, 
распределенной в зонах 3, 4, 5 по оси 7 (направленной к центру Земли) где М -— масса 
одной кубической ячейки. При этом в расчетах считается, что однородная масса с 
плотностью р, а аномальная - 1,2.-р. 
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агеа$ ш фе эгауцайопа| ехр]огайоп сотр!етеше4 Бу Фе гези!$ оЁ а загуеу шзйлитепе оЁ 
1$ апотаи$ 20пе сотр!етешагу ({есбпоо2?1ез (тазпейс уапайоп ехр]огайоп, 
зеосВеписа| ехр]огайоп затр!ез оЁ аш ап \умег ш Фе 50| сотропеп оЁ Ше апота!у 
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2опе). ТЬ$ агасе зитшаеа Бу Ше ргасйса! 1азК ш Йпатех тееогиез ог 511$ ОЕ уосашс 
аспуцу Бепеа 1 Фе 1се Беапп? сеорБузса! ап зеосВеплса| еЙес{$ ш зпо\у ог 1се 15 ИскКег, 
ап4 пеаг фе зи Расе оЁ Ше 1се соуег Потозепеой$ теЧпат (зисП аз зпо\ апА / ог 1се, аге 
Гурса! оЁ Ащагсйса). П Фете 1$ еу14епсе оЁ Ше отауцайопа| роепйа! аё Фе зи асе оЁ Фе 
Еацф сап Бе юип4 апота]ои$ 70пез, ипдег \у51сВ Феге аге ге?1оп$ \иф Чегет ШозрВегтс 
эгис@аге. ТВе сотршайопа! плоде! 15 а @гее-!ауег рагаПеерлре4 \Возе Базе агеа 1$ 20х20 
зпаге теегз ип! 5. Тре Бетор( оЁ еасЬ ]ауег 15 слуеп ш агбИгагу ипиз. ТБе изег Баз фе 
роз$Ииу ю сВоо$е Ше сотроз!йоп оЁ Ше 1ауег, ШегеБу свапеше Из депзиу, гезресиуеу. 

ш Фе патепйса]| ехрегитеп$ геуеае4 Фа Фе деулайоп$ оЁ Фе отауцанопа| Неа 
апотаНез дерепа оп фе 4ерй ап Фе гайо оЁ Из 4епзиу миф Фе 4епзиу оЁ Ще Тауег шт 
УВВ и ге$14е5. 


Надййила до редакций: 25.11.2016 


© В.Г. Писаренко, И.А. Варава, С.В. Корнеев, А.В. Кузько 143 


